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Abstract 

We extend to quenched QCD an earlier investigation of the chiral structure of the masses of 
the nucleon and the delta in lattice simulations of full QCD. Even after including the meson-loop 
self-energies which give rise to the leading and next-to-leading non-analytic behaviour (and hence 
the most rapid variation in the region of light quark mass), we find surprisingly little curvature 
in the quenched case. Replacing these meson-loop self-energies by the corresponding terms in full 
QCD yields a remarkable level of agreement with the results of the full QCD simulations. This 
comparison leads to a very good understanding of the origins of the mass splitting between these 
baryons. 



I. INTRODUCTION 



The quenched approximation is a widely used tool for studying non-perturbative QCD 
within numerical simulations of lattice gauge theory. With an appropriate choice of the 
lattice scale and at moderate to heavy quark masses, this approximation has been shown to 
give only small, systematic deviations from the results of full QCD with dynamical fermions. 
Although no formal connection has been established between full and quenched QCD, the 
similarity of the results has led to the belief that the effects of quenching are small and hence 
that quenched QCD provides a reasonable approximation to the full theory |]J . 

Improved lattice actions, together with advances in high performance computing, have 
been responsible for significant improvements in the calculation of baryon masses at moder- 
ate to light quark masses within the quenched approximation P, [| f|, ||. Simulations with 
dynamical fermions have proven to be more difficult, but results have been reported with 
pion masses as low as 320 MeV j|, |6| . 

The fact that one is restricted to quark masses much larger than the physical values 
means that, in addition to all the usual extrapolations (e.g., to the infinite volume and 
continuum limits), if one wants to compare with empirical hadron observables, one must also 
have a reliable method of extrapolation to the chiral limit. Any such extrapolation must 
incorporate the appropriate chiral corrections, arising from Goldstone boson loops, which 
give rise to rapid, non-linear variations as the chiral limit is approached. The importance of 
incorporating such behaviour has been successfully demonstrated for a number of hadronic 
observables, including masses |7|, [J, the sigma commutator ||, magnetic moments |TD|, [II], 



12, 1 3 1 , charge radii [14] and parton distribution functions [15, 16 



The impressive results found using these methods have led us to the present investigation 
of the problem of the chiral extrapolation of baryon masses in quenched QCD. The chiral 
properties within the quenched approximation are known to differ from those of full QCD 



in a number of very interesting ways [[17], [18], [19], |2(J, [H], p2| . For example, not only are 



the effective couplings at the pion-baryon vertices significantly altered in quenched QCD 
(QQCD) but, because the rj' behaves as a Goldstone boson in QQCD, one must also consider 
rj' loops. 

Here we first review previous work |7J which reported a successful method for extrap- 
olating baryon masses as calculated in full QCD lattice simulations. The modified chiral 
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structure of quenched baryon masses p0| is presented next. We show how to construct the 
various meson loop induced self-energies in order to preserve the leading-non-analytic (LNA) 
and next-to-leading non-analytic (NLNA) structure appropriate to QQCD, while incorpo- 
rating the established behaviour at heavier quark masses. This is followed by a detailed 
application to the extrapolation of the quenched N and A masses to the chiral limit. Fi- 
nally, we use the observed similarity of the structure of baryons stripped of their Goldstone 
boson clouds, in full and quenched QCD, to explore whether one can make a connection 
between the masses calculated in QQCD and those obtained in a dynamical simulation. 
The remarkable agreement obtained suggests a number of further tests and also leads us, 
with considerable confidence, to an interpretation of the origin of the N—A mass splitting. 

II. QCD EXTRAPOLATION 



In general, the coefficients of the LNA and NLNA terms in a chiral expansion of baryon 

(LI 
N 



masses are very large. For instance, the LNA term for the nucleon mass is 8m A) 



-5.6mJ (with and <5m^ NA ^ in GeV). With m, = 0.5 GeV, quite a low mass for current 
simulations, this yields 5m^ NA ^ = 0.7 GeV — a huge contribution. Furthermore, in this 
region hadron masses in both full and quenched lattice QCD are found to be essentially linear 
in m\ or equivalently quark mass, whereas 5m^ NA ^ is highly non-linear. The challenge is 
therefore to ensure the appropriate LNA and NLNA behaviour, with the correct coefficients, 
as m n — > 0, while making the transition to linear behaviour as increases, sufficiently 
rapidly to describe the actual lattice data. 

A reliable method for achieving all this was proposed by Leinweber et al. J7|. They fit 
the full (unquenched) lattice data with the form: 

M B = a B + (3 B ml + S B (m^,A), (1) 

where is the total contribution from those pion loops which give rise to the LNA and 
NLNA terms in the self-energy of the baryon. For the N these correspond to the processes 
N -> N-n -> N and N -> Avr -> N, while for the A we need A -> Att -> A and 
A — > Ntt — > A. Explicitly, 

7r | n 

N — ONN ~t a NAi 

= ^AA + a AN- ( 2 ) 



In the heavy baryon limit, these four contributions (B — > B'n — > S) can be summarised as: 

3 r r Hk w(fc) m 

where cj(/c) = a/& 2 + m 2 is the intermediate pion energy and curb' = {M B > — M B ) is the 
physical baryon mass splitting. The coefficients G BB * are standard SU(6) couplings and 
are summarised in Section 0. The ultraviolet regulator, u(k), has a very natural physical 
interpretation as the Fourier transform of the source of the pion field. The LNA and NLNA 
structure of these diagrams is associated with the infrared behaviour of the corresponding 
integrals and hence is independent of the choice of regular isat ion scheme. The use of a form 
factor effectively suppresses the self-energies like A 2 /m 2 for ^> A, the characteristic mass 
scale of the cutoff. A common choice of regulator, which we use throughout this work, is 
the dipole form, u{k) = A 4 /(A 2 + k 2 ) 2 . 

The linear term of Eq. (|J), which dominates for ^> A, models the quark mass depen- 
dence of the pion-cloud source — the baryon without its pion dressing. This term also serves 
to account for loop diagrams involving heavier mesons, which have much slower variation 
with quark mass. Given the current state of the art in lattice simulations, data in the low 
to intermediate mass range is unable to reliably constrain the parameter A. There is con- 
siderable phenomenological support for choosing a dipole mass parameter somewhat smaller 



than found for the axial form factor of the nucleon, which is 1.03 ± 0.04 GeV 23, El, 25 



However, it is important to understand that the anticipated development of supercomputing 
resources and techniques are such that A may be constrained by full QCD simulation data 
within the next five years. 

Fitting lattice results to Eq. (|l|) is straightforward. Upon calculating the described self- 
energies for a given choice of A, the fitting procedure amounts to a simple linear fit in 
and (3 B . 



III. QUENCHED CHIRAL PERTURBATION THEORY 

Standard chiral perturbation theory (xPT) is a low energy effective field theory built 
upon the symmetries of QCD [EJ], E7|| . It amounts to an expansion of Green's functions in 
powers of momenta and quark mass about the chiral limit (m q = 0). In the case of baryon 
masses, xPT tells us the leading behaviour of the quark mass expansion. Because xPT is 
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TABLE I: Coefficients of the lowest order non-analytic terms in the chiral expansions of the N 
and A masses, with values from both full and quenched QCD listed for comparison {v^ 1 = 167r/ 2 , 
AM = M A - M^). 

an effective field theory, the renormalisation procedure must be performed order by order 
in perturbation theory. At higher and higher order, more and more unknown parameters 
are introduced. These parameters only play a role in analytic terms of the expansion. The 
coefficients of the leading non-analytic terms are constrained by chiral symmetry |28|] — they 
are independent of regular isat ion and the order of the chiral expansion. In connecting the 
results of lattice QCD to the physical world it is essential that one incorporates the correct 
non-analytic structure of the low energy theory. 

Quenched (Q%PT) provides the analogous low energy effective theory for QQCD 



|T8| , [19], g0|. Sea quark loops are removed from QCD by including a set of degenerate, 
commuting (bosonic) quark fields. These bosonic fields have the effect of exactly cancelling 
the fermion determinant in the functional integration over the quark fields. This gives a 
Lagrangian field theory which is equivalent to the quenched approximation simulated on the 
lattice. The low energy effective theory is then constructed on the basis of the symmetries 
of this Lagrangian. 

The leading chiral expansion of baryon masses in the quenched approximation has been 
calculated by Labrenz and Sharpe |[20|| . For the reasons already mentioned in the Intro- 
duction, it differs from the corresponding expansion in full QCD. In particular, the chiral 
expansion coefficients take different values and new, non-analytic behaviour is also intro- 
duced. The explicit form can be expressed as 
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(a) (b) 

FIG. 1: Quark flow diagrams for the chiral rf loop contributions appearing in QQCD: (a) single 
hairpin, (b) double hairpin. 

M B = il4 0) + cf m n + c%m 2 n + cf ml + cf m£ + c B AL ml log m n + . . . (4) 
with the coefficients of the terms which are non-analytic in the quark mass listed in Table |. 



We note that in Ref. [20] the N and A were treated as degenerate states in the chiral 
limit. Experience in other situations suggests that it is more accurate to retain a finite mass 
difference, in which case off-diagonal terms such as N — > Air — > N lead to the non-analytic 
behaviour of the form m\ logm^. 

The contribution linear in m n is unique to the quenched approximation. The quenched 
theory therefore exhibits a more singular behaviour in the chiral limit. The origin of this 
behaviour is the Goldstone nature of the r/ in QQCD and specifically the process shown in 
Fig. |I](b). We note also that the coefficients of the chiral expansion involve new couplings, 
7 and 7', which are related to the flavour-singlet, hairpin-baryon couplings for iV and A 
respectively, illustrated in Fig. |l](a). In the formalism of Ref. |2(J these are related to the 
couplings of full QCD via the relations, 

J = D-F, 7' = 0, (5) 

as described in Appendix |B[ There is some uncertainty over the flavour singlet couplings, 
especially in connection with OZI violation associated with the U(l) axial anomaly p9 



While this may modify our calculated curves at extremely light quark mass, it would have 
no significant effect on the fit to lattice data at large quark mass nor on the comparison of 
current quenched and full QCD data. 
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TABLE II: Chiral couplings appearing in the self-energy integrals, Eq. (0) for full QCD and Eq. 



for QQCD. In numerical calculations we have used the couplings arising from SU(6) relations [30], 
C = —2D and H = —3D, with the tree-level parameters, F = 0.50 and D = 0.76. 



IV. QUENCHED SELF ENERGIES 

Our aim is to apply a similar procedure for the chiral extrapolation of quenched QCD 
data to that which has proven successful for the physical theory. That is, we wish to 
generalise Eq. ([!]) to replace the LNA and NLNA self-energy terms arising in full QCD by 
their quenched analogues. The pion loop contributions have the same kinematic structure 
as those in full QCD. A simple redefinition of the couplings, Gbb>, m the expressions for the 
self-energies ensures that the correct LNA and NLNA of QxPT is maintained. Thus, the 
analytic expressions for the pion cloud corrections to the masses of the N and A are of the 
same form as the full QCD integrals (c.f. Eq. (|3])): 

a * B ' - -m^l GBB ' L dk u(k)(u BB , + u(k)y (6) 

where the quenched couplings, Gbb', are listed in Table [H], together with their physical 
counterparts. 

Within the quenched approximation rf loops also contribute to the low energy effec- 
tive theory, whereas they are usually neglected in the physical case. This is because a 
re-summation of internal loop diagrams (coming from the fermion determinant) means that 
the rj' remains massive in the chiral limit of full QCD. On the other hand, the absence of 
these virtual loops in the quenched approximation causes the flavour singlet rj' to behave 
as a Goldstone boson [|l^, |19| . As a consequence of this feature of the quenched theory, 
there are two new types of loop contributions to be considered. A schematic view of these 
processes is shown in Fig. [l|. 

The first of these two contributions, shown in Fig. |I](a), arises from a single "hairpin" 
interaction. As discussed above, it is responsible for the term proportional to 7 (7') in the 
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TABLE III: Couplings used in flavour singlet rj' self-energies. We take m 2 , = 0.42 GeV , lying 



between phenomenological and lattice estimates [31, p2, 33 [. The momentum dependence of the 



double hairpin vertex, which is believed to be small, is neglected. 

chiral expansion of the iV (A) mass. These couplings are discussed in considerable detail in 
Appendix ||. The structure of this diagram is exactly the same as the pion loop contribution 
where the internal baryon is degenerate with the external state. The integral representing 
this diagram is then the same as that for cr BB , 

^>._ 3 N ^r dk ^m. (7) 

16vr 2 / 2 Jo u 2 (k) 
The factors N B \ providing the correct non-analytic behaviour in the chiral expansion 
(Eq. ([|)), are displayed in Table [TO . 

The second of these new r/ loop diagrams arises from the double hairpin vertex, pictured 
in Fig. [l|(b). This contribution is particularly interesting because there are two meson 
propagators and it is therefore responsible for the non-analytic term linear in — this 
term being unique to the quenched case. The integral corresponding to this self energy can 
be written in a similar way: 

^ ,(2) - 3 rdt k ' u2{k) (8) 



(2) 

Note the sign change and the higher power of oo in the denominator. The coefficients, N B , 
providing the correct non-analytic behaviour in Eq. (|J) — in this case the coefficient of 
— are given in Table |TTT| . 

The sum of these four contributions then gives the net meson-loop induced self-energies 
within the quenched approximation, 

- cr BB + cr BB , + a B +cr B . (9) 

The individual contributions to the A" and A masses over a range of pion mass are 
plotted in Figs. ^| and |3[ These are all evaluated with the dipole regulator mass parameter 
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FIG. 2: Various self-energy contributions to Mjy for dipole mass, A = 0.8 GeV. From top down at 



0.1 GeV 2 , the curves correspond to (where a ~ over the symbol denotes a quenched QCD 



n' (2) n' (1) 

contribution) , aj^ , o"^^, o^A' total quenched Sjy, crJj A , crJ^ N and total physical Ejy- 



A = 0.8 GeV. The corresponding self-energies from full QCD are also shown for comparison. 
We note that in QQCD the contributions are typically quite a bit smaller and the double- 
hairpin graph, a^ 2 \ is repulsive. The differences are enhanced for the A where a\ N is also 
repulsive. We observe that the rapid, non-linear behaviour (which is effectively much larger 
in full QCD) is restricted to the region m 2 < 0.2 GeV 2 , above which the self-energies are 
quite smoothly varying functions of the quark mass. 



V. FITTING PROCEDURE 



The lattice data considered in this analysis comes from the recent paper of Bernard et al. 
Q. These simulations were performed using an improved Kogut-Susskind quark action, 
which is known to have good scaling properties f34|| . Unlike the standard Wilson fermion 
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FIG. 3: Various self-energy contributions to Ma for dipole mass, A = 0.8 GeV. From top down at 
ml = 0.1 GeV 2 , the curves correspond to (where a ~ over the symbol denotes a quenched QCD 
contribution) a\ N , <r}} 2 \ total quenched Ea, ^ajv"' ^AA' ff AA an< ^ total physical In- 
action, masses determined at finite lattice spacing are excellent estimates of the continuum 
limit results. 

We are particularly concerned with the chiral extrapolation of baryon masses and how 
their behaviour is affected by the quenched approximation. In such a study, it is essential 
that the method of scale determination is free from chiral contamination. One such method 
involves the static-quark potential. As low-lying pseudoscalar mesons made of light quarks 
exhibit negligible coupling to hadrons containing only heavy valence quarks, the low energy 
effective field theory plays no role in the determination of the scale for these systems. In 
fixing the scale through such a procedure one constrains all simulations, quenched, 2-flavour, 
3-flavour etc., to match phenomenological static-quark forces. Effectively, the short range 
(0.35 ~ 0.5 fm) interactions are matched across all simulations. 
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A commonly adopted method involving the static-quark potential is the Sommer scale 



35| , |36fl . This procedure defines the force, F(r), between heavy quarks at a particular length 
scale, namely r — 0.5 fm. Choosing a narrow window to study the potential avoids compli- 
cations arising in dynamical simulations where screening and ultimately string breaking is 
encountered at large separations. The lattice data analysed in this report uses a variant of 
this definition, choosing to define the force at t\ = 0.35 fm via r\F{ri) = 1.00 [ffl. 

As we remarked earlier, the non-analytic chiral behaviour is governed by the infrared 
regions of the self-energy integrals. The fact that the lattice calculations are performed on a 
finite volume grid means that the self-energy integrals implicit in current lattice simulations 
do not include the exact chiral behaviour. It is important to take this into account in the 
fitting procedure and we therefore follow Ref. || in replacing the continuum self-energy 
integrals used in the fitting process by a discrete sum over the meson momenta available on 
the lattice: 

47T jH k 2 dk = J d 3 k*^(^pj ■ ( 10 ) 

fey fez 

The self-energy integrals calculated in this way are what should be directly compared with 
the lattice data, and we illustrate these by open squares in subsequent figures. Upon obtain- 
ing the optimal fit parameters, one can evaluate the integrals exactly and therefore obtain 
the infinite- volume, continuum limit. The latter is the result which should be compared 
with experiment at the physical pion mass. 

We now proceed to fit quenched lattice data with the form 

M B = a B +(3 B ml + £ fl (m w ,A), (11) 

(by analogy with the form used in full QCD, Eq. ([I])), with the self-energies evaluated, as we 
have just outlined, using the momentum grid corresponding to the specific lattice simulation. 
The linear terms in Eq. (|II]) may be thought of as accounting for the quark mass dependence 
of the pion-cloud source. This form then automatically includes the expected heavy quark 
behaviour where the tt and rj' loop contributions are suppressed. 

The form factor, which models the physical structure of the meson-baryon vertex, char- 
acterises the finite size of the pion source. Quenched simulations of hadronic charge radii 



performed at moderate to heavy quark masses [37] have been demonstrated to be consis- 



tent with experiment once the meson-cloud properties of full QCD are taken into account 



14] , |38fl . This indicates that the size of the meson-cloud source is expected to be of similar 

11 



1.2 



1.0 



0.8 - 



<o 0.6 



0.4 - 



0.2 



0.0 




0.65 0.70 



0.75 0.80 0.85 
A (GeV) 



0.90 0.95 



FIG. 4: The value 5 is a measure of the difference between the quenched and dynamical data sets 
after accounting for the relevant self-energy diagrams. This measure is proportional to the net area 
contained between the straight lines obtained from the fits and has been normalised to the case 
where the self-energy diagrams are totally neglected. 



size in both quenched and physical QCD. For this reason we proceed to fit both quenched 
and physical data with a common value of A. 

For a fixed choice of A, fitting to lattice data amounts to a linear fit in a and (3. It 
turns out that, for a range of values of A, the values of a and (3 found for the QQCD 
data are surprisingly close to the values found for the fit to dynamical QCD data. This 
strongly suggests that the self-energies included here, which contain the LNA and NLNA 
behaviour appropriate to each type of simulation, contain the primary effect of quenching. 
To illustrate the point, Fig. |] shows a measure, 5, of the difference between the quenched 
and dynamical data sets over the range of m n considered. This measure is proportional 
to the net area contained between the straight lines obtained from the fits and has been 
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FIG. 5: Fit (open squares) to lattice data 0] (Quenched A, Dynamical A) with adjusted self- 
energy expressions accounting for finite volume and lattice spacing artifacts. The infinite-volume, 
continuum limit of quenched (dashed lines) and dynamical (solid lines) are shown. The lower 
curves and data points are for the nucleon and the upper ones for the A. 



normalised to the case where the self-energy diagrams are totally neglected. The improved 
agreement between data sets over the range of dipole masses highlights the effectiveness 
of this self-energy correction. It is also worth noting that the x 1 /dof is also improved by 
incorporating the self-energies into the fit. For the preferred dipole mass, A = 0.8 GeV, this 
is better by a factor 2. 

Results of both the physical and quenched fits are shown together in Fig. [5| The pa- 
rameters of the best fits are displayed in Table [IV]. Here we see the remarkable agreement 
of the linear term of our fitting formulae, Eqs. ([!]) and (pT|). This strongly suggests that 
the behaviour of the meson-cloud source is very similar in quenched and full QCD. The 
primary difference between the quenched and physical results can then be described by the 
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TABLE IV: Best fit parameters for both full and quenched data sets with dipole regulator, A = 
0.8 GeV. The second set correspond to a simple linear fit, where the self-energy contributions have 
been neglected. All masses are in GeV. 



meson-loop induced self-energies. 

This observation suggests that it may well be possible to make a connection between 
quenched simulations and hadron properties in the real world. One would fit quenched data 
with appropriate self-energies to obtain the linear behaviour of the meson-cloud source. 
Then the quenched self-energies would be replaced by their full-QCD counterparts, hence 
obtaining more physical results. It is clearly very important to test this result further on 
other hadrons (e.g. for other members of the octet) and against dynamical simulations at 
lower quark masses. 

VI. A-N HYPERFINE SPLITTING 

The analysis of lattice data has demonstrated the ability to describe the primary differ- 
ence between quenched and dynamical simulations in terms of the meson-loop self-energies. 
Figure ||] shows the difference in the self-energy terms for the N and A in quenched and 
full QCD, for several values of the common dipole mass. It is quite clear that there is a 
difference of between 150 and 250 MeV between the quenched and full QCD cases. Since 
this difference was essential in accounting for the clear differences in the behaviour of the 
baryon masses in QQCD and full QCD shown in Fig. |], we have some confidence in using 
these results to say how much of the physical N-A mass splitting is associated with pion 
loops and how much comes from short range processes, such as gluon exchange. In fact, an 
examination of Fig. ^|for the case of full QCD suggests fairly clearly that only about 50 MeV 
of the observed 300 MeV N-A splitting arises from pion loops. Of course, this result is more 
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FIG. 6: Meson-loop contribution to the A-N mass splitting in both quenched and full QCD - for 
several values of the dipole mass. 



dependent on the assumption of the same dipole mass parameter at every vertex than the 
fits to the N and A masses individually. Nevertheless, it seems unlikely that more than a 
third of the total splitting could come from this source. 



VII. CONCLUSIONS 



We have investigated the quark mass dependence of the N and A masses within the 
quenched approximation. The leading chiral behaviour of hadron masses is known to differ 
in quenched QCD from the physical theory. This knowledge has been used to guide us in 
the construction of a functional form which both reproduces the correct chiral structure, 
and is consistent with current lattice simulations. This procedure of fitting lattice data 
with a linear term together with the meson loop corrections which give rise to the LNA 
and NLNA behaviour has been demonstrated previously to fit dynamical QCD simulation 
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results remarkably well. Here we have shown that the application of the same procedure 
to quenched results is able to consistently fit the data in that case as well. Remarkably, a 
comparison of the two fits suggests that the properties of the N and A, stripped of their pion 
clouds, are essentially the same in quenched and full QCD, once the scale is set using the 
Sommer scale appropriate to heavy quark systems. It is clearly essential to test this finding 
against further full QCD simulations at lighter quark masses as well as for other hadrons. 

We have demonstrated that although the quenched approximation gives rise to more 
singular behaviour in the chiral limit, this is not likely to be observed in lattice simulations 
as these contributions are quickly suppressed with increasing quark mass. Indeed our results 
suggest that it will be very hard to detect any significant chiral curvature in the case of the 
nucleon, while for the A there may be some small, upward curvature. The A-iV mass 
splitting increases to around 400 MeV at the physical point in QQCD. As a consequence of 
this behaviour, the A mass in the quenched approximation is expected to differ from the 
physical mass by approximately 25%. Finally, we have shown that while a fraction of the 
physical iV-A mass difference can be attributed to a difference in pion self-energy loops, 
this is unlikely to amount to more than a third of the observed splitting. 
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APPENDIX A: ANALYTIC INTEGRATION 

Here we summarise analytic expressions for the self-energy integrals. It should be noted 
that these expressions are not used in fitting lattice data. For the purpose of fitting, the 
continuum integral is replaced by a discrete sum over the available momenta on the corre- 
sponding lattice, as described in Eq. ([10|). 

Firstly we consider the case of the simple meson-loop digram where the internal baryon 
line has degenerate mass with the external state, 




(Al) 
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Using a sharp cut-off, u(k) = 9 (A — k), the integral can be expressed as 

f , / A \ A 3 A A 

arctan — M Am^ . (A2) 



16tt 2 ^ V w \m. 
The LNA behaviour of this can then be immediately read from this, 



I LNA 3G 3 
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Alternatively, our preferred dipole u(k) = A 4 / (A 2 + /c 2 ) 2 , also provides an analytic expression 
for this self-energy, 

3G A 5 



a 



(ml + AAm^ + A 2 ) . (A4) 



512tt/ 2 (A + m 7r ) 4 

This gives precisely the same LNA behaviour as the sharp cut-off, as expected because the 
non-analytic behaviour is due to the infrared behaviour of the integral. It is associated with 
the residue of the pion propagator pole, and hence independent of an ultraviolet cut-off. 
The case of the double meson propagator can also be performed analytically, 

3 Ar f°° J7 k 4 u 2 (k) 

For u(k) =9(A-k), 



q 3774A + 2A 3 - ZmJml + A 2 ) arctan (- 

° = T7TTT2 N Wr^~^ — " (A6) 

16tt 2 / 2 2(m 2 +A 2 ) 



For u(k) = A 4 /(A 2 + A; 2 ) 2 



_ 3N A 5 (m^ + 5A) 

° 5127T/ 2 (A + 777^ ■ 1 } 



Once again both integrals give the same LNA behaviour, 

I LNA 9^ /* x 

For the off-diagonal contributions, where the internal baryon is not degenerate with the 
external state, 



; G I dk ( A9 ) 



167T 2 / 2 J u(k)(u B B> +uj(k)) : 
with ujbb 1 finite. The results do not have a simple form. The full expression for the case of 
a sharp cut-off form factor can be found in Ref. j7|. We show the LNA contribution to this 
diagram for reference, 



lna 9G mi 



a\^ A = ~ * log m^. (A10) 
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APPENDIX B: FLAVOUR SINGLET IN FULL AND QUENCHED QCD 



This appendix serves to clarify the derivation of the hairpin meson-baryon couplings in 
quenched xPT. 

The flavour singlet rj' remains light in the quenched approximation, and is therefore an 
effective degree of freedom in the low energy sector. Such excitations must therefore be 
incorporated into the low-energy analysis. Within full QCD, resummation of internal loop 
diagrams renders the rf massive and hence it plays no role in the low-energy dynamics. For 
this reason couplings to such flavour singlet states are neglected. In our analysis, we wish 
to compare the low-energy structure of the quenched and physical theories. In this case, a 
flavour singlet coupling, like NNrj', must be included in the chiral Lagrangian of full QCD 
in order that it is treated on equal footing with the quenched theory. This coupling will not 
alter any results of the physical theory as any diagram would involve the propagation of a 
heavy rf . 

Here we derive best estimates for the flavour singlet couplings in quenched QCD. This 
is achieved by comparison of the quenched and full chiral Lagrangians, under the standard 



assumption that the couplings exhibit negligible change between the two theories PU|. We 



follow the notation of Labrenz and Sharpe p(J in the analysis of such contributions. All 



symbols retain the same meaning, unless otherwise specified. 
The chiral Lagrangian for full QCD can be expressed as 

£ = Cir + Cbtt + £tit ■ (Bl) 

The standard octet and decuplet Lagrangians are given with an additional coupling to an 
SU(3) flavour singlet state. For clarity we label the octet and singlet parts of the meson 
matrix, A: 

C B n = itr(Bv ■ VB) 

+2£>tr (BS^{A°f, B}) + 2Ftr {BS^Af, B}) 

+2/jb D tr (B{M + , B}) + 2/iMr (B[M + , B]) 

+2/A,tr (BB)ti (M + ) + 2g s tr (BS^B)tx (A s j n ) , (B2) 

Cm = —iT u (v ■ V)T V + AMT U T U 

+2Hf v S ti AfT v + C{f v AfB + BA° v ct T u ) 
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+cT u M + T u - aT u T v tr (M + ) 

+2g' s T v S»T u ti {A s p . (B3) 

The new parameters, g s and g' s , describe couplings of the flavour singlet 7/ to baryon octet 
and decuplet states respectively. Within full QCD the single vertex has two topologically 
different quark flow diagrams as illustrated by the left and right-hand vertices of Fig. |l|(a). 
The left is that of a qq insertion on one of the valence quark lines and the right is a pure 
gluonic coupling through a hairpin-style qq annihilation. The total coupling is a sum of 
these two contributions. Denoting the hairpin vertex coupling by 7qcd and 7q CD for octet 
and decuplet baryons respectively we have 

9s = -j= 9 V 'NN + 7qcd , (B4) 

9's = 4= + 7QCD ■ (B5) 

The first of these interactions, g v >NN (fly a a) is related to the axial couplings by SU(6) 
phenomenology. We take the standard approach and assign 

Qrj'NN = V2g vNN = \J\{3F - D) , (B6) 
fly A A = V2g vAA = • (B7) 

The effective chiral Lagrangian of quenched QCD is [^Oj 

CM = 4 Q) + Cfi + , (B8) 

where meson and baryon states are now understood to be constructed of ordinary quarks 
and bosonic quarks. The general Lagrangian for the heavy fields can be written in terms of 



the rank-three tensor fields as defined in Ref. 20 , B and T 



4J = i(Bv-VB) 

+2a{BS^ l BA fl ) + 2[3(BS tl A IJt B) 
+2 ls {BS»B) str(A M ) + a M (BBM + ) 

+Pm(BM + B) + a(BB) str(7W+) , (B9) 
Cfl = -i(T u (v ■ V)%) + AM(f u %) 



+2H(T V S"A I X) - ^C[f"A u B + BA V T 
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+2i s (T v S»%) str(^) 

+cT v M + % - a(T v T v )ti (M + ) . (BIO) 

It should be noted that the terms 7 S and 7^ describe both types of flavour-singlet coupling, 
not just that arising through the hairpin alone. Similarly to Eq. ( P5|) , in the quenched 
theory these can be described by 

1 

Is = -j=g v >NN + l, (Bll) 

7; = i^AA + V, (B12) 

where the terms 7 and 7' now correspond to the pure hairpin couplings as used in Ref. |2D| . 
Here we also note the the terms g v 'NN and g v >AA are unchanged in going to the quenched 
theory, this is consistent with the assumption that the chiral parameters F and D are 
unchanged between the two theories. 

One can then relate the quenched chiral Lagrangian back to that of full QCD by restricting 
the indices on the tensor fields, B and T, to those corresponding to the physical quarks. 
The details of this procedure are described in Ref. p0| . Performing these restrictions on the 
octet-baryon, quenched chiral Lagrangian (Eq. |B9| ) one finds: 

4*1* = i{Bv-VB)\ R 

+ |(2« - /3)tr (BS^A^B) + i(-a - 4/3)tr (BS^BA^) 

+~(a + 4(3 + 6 7s )tr (BS^B)tr (AJ 

+~(2a M - P M )ti (BM+B) + U-a M - Af3 M )tx (BBM+) 

6 O 

+\{a M + A(3 M + 6a)tr (BB)tx (M + ) . (B13) 


Equating this with Eq. (|B!2|) gives: 



|(2a-0) 


= 2D + 2F, 


(B14) 




= 2D-2F, 


(B15) 


l -(a + A[3 + Q ls ) 


= 2g s , 


(B16) 


^(2a M - 13m) 


= 2fib D + 2\ib F , 


(B17) 


-{-a M - 4:(3 M ) 



= 2\ibr) — , 


(B18) 


-(a M + 4f3 M + 6cr) 


= 2{ib . 


(B19) 
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In extracting the flavour-singlet part, Eq. ( |B16|) provides us with 

-a + -{3 + J - g v > NN + 27 = J - g v , NN + 27 QCD , (B20) 



and combining with Eqs. (B14.B15) one arrives at 

j = 1qcd + D-F. (B21) 
The restrictions are much simpler for the decuplet case and one finds 

l' = 7qcd • (B22) 

In estimating the hairpin-type couplings in full QCD one assumes that they are relatively 
small, 7qcd *C g v 'NN, due to OZI-type suppression |39|]. With analogous arguments for the 
decuplet, we take 7qcd — 7qcd = 0- We do n °te that the U(l) axial anomaly may be 
effective in overcoming the OZI rule in the case of 77' couplings |29|], but as we mentioned 
in the text the main conclusions of our present analysis are not very sensitive to the precise 
value of the r/'-nucleon coupling. 
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